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Abstract 

A uniform extracellular stimulus triggers cell-specific patterns of Ca 2+ signals, even in genetically identical cell populations. 
However, the underlying mechanism that generates the cell-to-cell variability remains unknown. We monitored cytosolic 
inositol 1,4,5-trisphosphate (IP 3 ) concentration changes using a fluorescent IP 3 sensor in single HeLa cells showing different 
patterns of histamine-induced Ca 2+ oscillations in terms of the time constant of Ca 2+ spike amplitude decay and the Ca 2+ 
oscillation frequency. HeLa cells stimulated with histamine exhibited a considerable variation in the temporal pattern of 
Ca 2+ signals and we found that there were cell-specific IP 3 dynamics depending on the patterns of Ca 2+ signals. RT-PCR and 
western blot analyses showed that phospholipase C (PLQ-fil, -(33, -(54, -yl, -83 and -e were expressed at relatively high levels 
in HeLa cells. Small interfering RNA-mediated silencing of PLC isozymes revealed that PLC-(31 and PLC-f54 were specifically 
involved in the histamine-induced IP 3 increases in HeLa cells. Modulation of IP 3 dynamics by knockdown or overexpression 
of the isozymes PLC-pM and PLC-fS4 resulted in specific changes in the characteristics of Ca 2+ oscillations, such as the time 
constant of the temporal changes in the Ca 2+ spike amplitude and the Ca 2+ oscillation frequency, within the range of the 
cell-to-cell variability found in wild-type cell populations. These findings indicate that the heterogeneity in the process of IP 3 
production, rather than IP 3 -induced Ca 2+ release, can cause cell-to-cell variability in the patterns of Ca 2+ signals and that 
PLC-fSI and PLC-(54 contribute to generate cell-specific Ca 2+ signals evoked by G protein-coupled receptor stimulation. 

Citation: Ishida S, Matsu-ura T, Fukami K, Michikawa T, Mikoshiba K (2014) Phospholipase C-pi and P4 Contribute to Non-Genetic Cell-to-Cell Variability in 
Histamine-induced Calcium Signals in HeLa Cells. PLoS ONE 9(1): e86410. doi:10.1371/journal.pone.0086410 

Editor: Laszlo Csernoch, University of Debrecen, Hungary 

Received August 9, 2013; Accepted December 12, 2013; Published January 27, 2014 

Copyright: © 2014 Ishida et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits 
unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. 

Funding: This work was supported by grants from the Ministry of Education, Culture, Sports, Science and Technology of Japan to T. Matsu-ura (22770227), T. 
Michikawa (20370054), and K. Mikoshiba (2022007). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of 
the manuscript. 

Competing Interests: The authors have declared that no competing interests exist. 
* E-mail: t-michikawa@brain.riken.jp (TM); mikosiba@brain.riken.jp (KM) 

na Current address: Laboratory for Cell Signaling Dynamics, RIKEN Quantitative Biology Center, Osaka, Japan 

nb Current address: Department of Molecular and Cellular Physiology, University of Cincinnati, Cincinnati, Ohio, United States of America 



Introduction 

Many extracellular stimuli trigger increases in the cytosolic 
concentration of Ca 2+ ([Ca 2+ ]) that control a wide range of 
physiological processes, including fertilization, proliferation, de- 
velopment, learning and memory, contraction, and secretion [1,2]. 
In a multitude of cell types, released Ca 2+ causes oscillatory 
changes in [Ca 2+ ], and the frequency of Ca 2+ oscillations is 
correlated with the stimulus intensity [3-5] , while the time course 
of an individual Ca 2+ spike depends on the type of receptor 
stimulated [6,7]. The versatility of Ca 2+ signals is considered to be 
dependent on the spatiotemporal patterns of intracellular Ca + 
signals [8]. 

Extracellular stimuli, such as growth factors, hormones, and 
neurotransmitters, activate an intracellular effector molecule, 
phosphoinositide-specific PLC, that catalyzes the conversion of 
phosphatidylinositol 4,5-bisphosphate (PIP2) to inositol 1,4,5- 
trisphosphate (IP3) and diacylglycerol [9]. IP3 acts as a second 
messenger in many cell types, where its primary effect is to trigger 



the release of Ca + ions from intracellular stores and generate 
cytosolic Ca + signals. The forms of the Ca + signals vary from cell 
to cell, even in clonal populations of genetically identical cells, 
while individual cells show characteristic and reproducible Ca + 
responses in terms of amplitude, frequency, and shape of the 
transient as well as the latency of onset of the initial Ca 2+ rise [5]. 
Therefore, Ca + signals have been referred to as cell "finger- 
prints". 

Real-time monitoring of cytosolic IP 3 in living cells is pivotal for 
understanding the mechanism that underlies the generation of cell- 
specific patterns of Ca 2+ signals. In a previous study, we monitored 
the cytosolic IP 3 concentration ([IP3]) changes during agonist- 
evoked Ca 2+ oscillations in HeLa cells using a genetically-encoded 
fluorescent IP 3 sensor IP 3 R-based IP 3 sensor 1 (IRIS-1) [10]. The 
observed IP 3 dynamics were unexpectedly complex because [IP 3 ] 
did not return to its basal level during the intervals between Ca + 
spikes, and IP 3 gradually accumulated in the cytosol with little or 
no fluctuation during cytosolic Ca 2+ oscillations [10]. Manipula- 
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tion of the cytosolic Ca + increases during agonist application 
indicated the presence of both Ca 2+ -dependent and Ca 2+ - 
independent components for IP 3 generation in this cell type 
[10]. The molecular natures of these components and, more 
fundamentally, the roles of IP3 dynamics in the generation of the 
cell-specific patterns of Ca 2+ signals remain elusive. 

There are six families of mammalian PLC enzymes (PLC-P, -y, 
-5, -s, -'(,, and -1"|) consisting of 13 isoforms in humans [11,12]. 
PLC-P is activated by the Got and Gfiy subunits of heterotrimeric 
G proteins downstream of G protein-coupled receptors (GPCRs). 
PLC-y is activated through receptor and non-receptor tyrosine 
kinases. PLC-8 is activated by cytosolic Ca + and the high- 
molecular-weight G protein Gh. PLC-s is regulated directiy by 
small GTPases from the Ras family, and its enzymatic activity can 
be stimulated by subunits of heterotrimeric G proteins and small 
GTPases from the Rho family. PLC-i^ was identified as a sperm- 
specific PLC, and shows extremely high Ca 2+ sensitivity for its 
enzymatic activity compared with the other PLC isoforms. PLC-r) 
is suggested to be activated through GPCR stimulation. Although 
various tissues express multiple PLC isozymes, the effects of the 
individual PLC isozymes on the generation of Ca 2+ signals are 
poorly understood. 

In this study, we measured the IP3 dynamics using IRIS-1 in 
cells showing different patterns of Ca 2+ signals within a clonal 
population of genetically identical HeLa cells stimulated with 
histamine. The IP3 dynamics were characteristically different 
depending on the patterns of cellular Ca 2+ signals, suggesting that 
the heterogeneity in IP3 production contributes to the cell-to-cell 
variability observed in Ca 2+ signals. We combined an RNAi 
technique with the real-time IP3 imaging technique to identify the 
PLC isozymes involved in the generation of IP :i production in this 
cell type, and found that at least two isozymes in the PLC-P family, 
PLC-Pl and PLC-P4, are involved in histamine-induced IP3 
increases in HeLa cells. RNAi-mediated silencing and overexpres- 
sion experiments for PLC-P 1 and PLC-P4 revealed that these 
isozymes differentially contribute to the determination of the 
characteristic parameters of Ca 2+ signals, i.e. the Ca 2+ oscillation 
frequency and the time constant of Ca 2+ spike amplitude decay. 
These findings provide important clues toward understanding the 
molecular mechanism that underlies the generation of cell-specific 
patterns of Ca 2+ signals in genetically identical cell populations. 

Results 

Reproducible Cell-specific Patterns of [Ca 2+ ] and [IP 3 ] 
Changes in HeLa Cells Stimulated with 3 j_i/Vl Histamine 

To examine whether the cell-specific patterns of Ca 2+ oscilla- 
tions were accompanied by characteristic IP 3 dynamics in 
individual cells, we performed simultaneous monitoring of 
[Ca 2+ ] and [IP3] at the single-cell level. Indo-5F was loaded into 
IRIS-1 -expressing HeLa cells and the cells were stimulated twice 
with 3 U.M histamine with an interval of approximately 20 min. 
Consistent with previous studies [5,13-15], there was marked 
intercellular heterogeneity in the Ca 2+ dynamics, which varied in 
the pattern (sustained oscillations or damped oscillations) and the 
frequency of Ca 2+ oscillations (for more details, see below), even 
within the same microscopic field of view (Fig. 1A). These cell- 
specific Ca 2+ patterns were reproducible upon repetitive histamine 
additions (Fig. IB), as observed in different cell types [5,13-15], 
suggesting that individual HeLa cells have characteristic Ca 2+ 
dynamics. As shown in Fig. IB, the patterns of IP3 dynamics also 
differed between cells. In cells showing sustained Ca 2+ oscillations 
(e.g. cell 1 in Fig. IB), IP 3 was maintained at relatively low levels 
during the period of stimulation, and fluctuations in [IP3] that 



synchronized with Ca + spikes were observed with almost constant 
amplitudes. In cells showing damped Ca 2+ oscillations (e.g. cell 2 in 
Fig. IB), the initial increase in [IP 3 ] was relatively high and the IP :i 
level gradually increased during the period of stimulation. [IP3] 
fluctuations that synchronized with Ca 2+ spikes were also observed 
in cells showing damped Ca + oscillations, but the amplitude of 
these [IP3] fluctuations gradually decreased depending on the 
amplitude of the Ca 2+ spikes. In both cell types, the IP 3 dynamics 
were reproducible (Fig. IB). These findings indicate that the IP3 
dynamics are characteristic of individual cells, and that the IP 3 
dynamics may be involved in the determination of cell-specific 
Ca 2+ responses. 

[Ca 2+ ] and [IP 3 ] Dynamics Observed in Cells Stimulated 
with Different Concentrations of Histamine 

As shown in Fig. 1, there were two prominent types of cells 
depending on their temporal patterns of [Ca 2+ ] changes after 
stimulation with 3 uM histamine. To investigate whether these 
differences arose from heterogeneity in the histamine sensitivity of 
the cells, as proposed by Rooney et al. [15], or from other 
components in the signal transduction pathway, we measured the 
[Ca 2+ ] and [IP 3 ] chan ges in the same cells stimulated with three 
different concentrations of histamine. Representative results are 
shown in Fig. 2A (cells showing sustained Ca 2+ oscillations at 
3 |XM histamine) and Fig. 2B (cells showing damped Ca 2+ 
oscillations at 3 uM histamine). To analyze the results quantita- 



A 




1 min 



Figure 1. Reproducible cell-specific patterns in Ca and IP 3 
responses in HeLa cells stimulated with histamine. (A) Venus 
fluorescence image of IRIS-1 -expressing HeLa cells before histamine 
stimulation. Bar, 20 am. (B) Changes in 460-510 nm emission of lndo-5F 
signals (F/F a where F 0 is the basal level of F; top) and ECFP/Venus 
emission ratio of IRIS-1 signals (AR/R n ; AR was defined as R - R 0 , where 
R 0 is the basal level of R; bottom) after repeated additions of 3 |oJVI 
histamine with a 20-min interval in two different cells within the same 
field of view shown in (A). The horizontal broken lines indicate the 
baseline levels of the IRIS-1 and lndo-5F signals. The vertical broken 
lines indicate the onsets of stimulation. 
doi:10.1371/journal.pone.0086410.g001 
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tively, we measured the time constants for the progressive decrease 
in the amplitude of Ca 2+ spikes (Fig. SI A). Sustained oscillations 
had larger time constants of Ca 2+ spike amplitude decay, while 
damped oscillations had smaller time constants of Ca 2+ spike 
amplitude decay. Fig. 2C shows the relationships between the time 
constants and the histamine concentrations applied. There were 
large cell-to-cell varieties for all three histamine concentrations 
examined, and the inverse time constants varied from less than 
0.001 s -1 (stable sustained oscillations) to 0.3 s _1 (rapidly damped 
oscillations). There was no obvious dependence of the time 
constant on the histamine concentration, and individual cells 
showed similar time constants irrespective of the histamine 
concentration (Fig. 2C). These findings suggest that a difference 
in histamine sensitivity is not the primary mechanism for the 
generation of different patterns of Ca 2+ signals and that the time 
constants are characteristic of individual cells. Because there were 
two peaks in the histogram of the time constants observed in cells 
stimulated with 3 uM histamine (Fig. SIB), we classified the cells 
into two groups, namely S-cells (cells showing sustained oscilla- 
tions: red in Fig. 2C) and D-cells (cells showing damped 
oscillations: blue in Fig. 2C), depending on the time constants of 
Ca 2+ amplitude decay observed at 3 uJVI histamine with a 
threshold of 0.0056 s . 

Fig. 2D shows the relationships between the histamine 
concentrations and the Ca 2+ oscillation frequencies in S-cells 
(red) and D-cells (blue). The frequency increased depending on the 
concentration of histamine, with slopes of 2.2±1.1 mHz/uM 
histamine (mean ± SD, n = 14) for S-cells and 1.2 ±1.0 mHz/uM 
histamine (mean ± SD, n = 23) for D-cells. There were large cell- 
to-cell variabilities in the Ca 2+ oscillation frequencies, but the 
frequencies in S-cells tended to be smaller than those in D-cells 
(Fig. 2D-G). We detected significant differences (P<0.05, 
Student's i-test) in the Ca 2+ oscillation frequencies between S-cells 
and D-cells following stimulation with 1 u,M histamine (mean ± 
SD: 20.5±10.3 mHz, «=11, for S-cells vs. 31.1 ± 10.6 mHz, 
n = 1 1, for D-cells) and 3 |iM histamine (mean ± SD: 
25.0±6.9mHz, n= 13, for S-cells vs. 32.5± 1 1.2 mHz, n = 23, 
for D-cells) (Fig. 2E and F). These findings indicate the presence of 
a correlation between the interval of Ca 2+ spikes and the time 
dependence of the change in the amplitude of Ca 2+ spikes. 

Representative traces showed that the dependence of the Ca 2+ 
spike amplitude on the histamine concentration differed between 
S-cells and D-cells (Fig. 2A and B). We quantified the peak 
amplitudes of the first Ca 2+ spike following stimulation with three 
different concentrations of histamine (Fig. S2A and B). The 
amplitude of the first spike of S-cells remained almost constant 
irrespective of the histamine concentration applied, while that of 
D-cells was dependent on the histamine concentration (Fig. S2B). 
These results suggest that the histamine concentration-dependence 
of the amplitude of the first Ca 2+ spikes is tighdy coupled to the 
rate of changes in the amplitude of the following Ca 2+ spikes in 
HeLa cells. 

Fig. 2H^J shows histograms of the integrals of the IP :i signals 
during stimulation with the three different concentrations of 
histamine. As shown in Fig. S2C, the mean values of the integrated 
IP3 signals were dependent on the histamine concentration in both 
S-cells and D-cells. The IP3 productions in S-cells were 
significandy smaller than those in D-cells under all three 
conditions, suggesting that the amount of IP3 production and/or 
IP3 dynamics should be involved in the formation of the cell- 
specific patterns of [Ca 2+ ] changes. To understand the relationship 
between IP3 dynamics and Ca 2+ dynamics more precisely, we 
investigated the changes in Ca 2+ dynamics by modifying IP 3 
production in the following analyses. 



PLC Isozymes Expressed in HeLa Cells and their Specific 
Knockdown Using Small Interfering RNAs (siRNAs) 

In a previous study [10], we found that there are multiple 
components of IP3 production in HeLa cells stimulated with 
histamine (see below). To identify the PLC isozymes involved in 
the generation of IP3 dynamics in histamine-stimulated HeLa cells, 
total RNA was isolated from HeLa cells and subjected to RT-PCR 
analyses with specific primer pairs (Table 1). We found that PLC- 
pT, -P3, -(34, -yl, -83 and -s were expressed at relatively high levels 
among the 13 PLC genes (Fig. 3 A). All of the isozymes were 
successfully detected by western blotting with isozyme-specific 
antibodies (Fig. 3B). In the western blot analyses, two spliced 
variants of PLC-fJl and PLC-P4 were detected (Fig. 3B). Next, we 
examined the effects of specific knockdown of the PLC isozymes 
using siRNAs. As shown in Fig. 3C, the expressions of the 
endogenous PLC isozymes, except for PLC-s, were markedly 
reduced by two different siRNAs for each isozyme. We confirmed 
that the PLC isozyme-specific knockdown did not affect the 
expression levels of the other PLC isozymes in HeLa cells (Fig. S3). 

Identification of PLC Isozymes Involved In Histamine- 
induced [IP 3 ] Increases in HeLa Cells 

There are multiple components in the IP3 increases observed in 
HeLa cells [10]: (1) a fast IP3 increase in response to histamine 
receptor stimulation without cytosolic Ca 2+ increase (Fig. 4A, left); 
(2) a slow IP3 increase in response to cytosolic Ca 2+ increase alone 
(Fig. 4A, middle); and (3) a complex [IP3] change composed with a 
fast transient (1st component) and a slow gradual increase (2nd 
component) in response to the combination of histamine receptor 
stimulation and cytosolic Ca 2+ increase (Fig. 4A, right). The 
presence of cytosolic Ca 2+ enhanced the rate of the histamine- 
induced IP 3 increase and also elicited a delayed negative feedback 
(Fig. 4B). We attempted to identify the PLC isozymes contributing 
to each component by measuring the IP3 dynamics in PLC 
isozyme-specific knockdown cells. To exclude the effects of Ca 2+ 
released through IP 3 receptor stimulation on the [IP3] changes, the 
cells were pretreated with thapsigargin to deplete intracellular 
Ca 2+ stores, and cytosolic [Ca 2+ ] was increased using the 
capacitive Ca 2+ entry pathway activated by Ca 2+ store depletion, 
using previously described methods [10]. We used the maximal 
IRIS-1 signal change or the increasing rate of the IRIS-1 signal as 
an index of PLC activity in the following experiments because 
these values were not affected by the IRIS-1 expression level, the 
resting level of [IP3], or the increasing rate of [Ca 2+ ] or total 
amount of [Ca 2+ ] increase mediated by capacitive Ca + entry (Fig. 
S4). 

Figure 4C shows the results of PLC isozyme knockdown on the 
[IP3] increase in response to histamine receptor stimulation alone. 
Among the five isozymes examined, only PLC-P4 knockdown 
caused significant impairments of the [IP3] increase irrespective of 
the siRNA used. The maximal IRIS-1 change was significandy 
decreased by 46% (siRNA#l) and 71% (siRNA#2) compared 
with the mean change with the control siRNAs mGC and hGC 
(Fig. 4C). These findings indicate that PLC-P4 is selectively 
activated downstream of histamine receptors in the absence of 
[Ca 2+ ] elevation in HeLa cells. 

Figure 4D shows the results of PLC isozyme knockdown on the 
[IP3] increase in response to cytosolic Ca 2+ increase alone, which 
was elicited by the application of extracellular 2 mM CaCl 2 
whithout histamine. In thapsigargin-treated HeLa cells, store- 
operated Ca 2+ influx channels on the plasma membrane are 
activated by the depletion of intracellular Ca 2+ stores and the 
addition of extracellular Ca 2+ immediately elicites [Ca 2+ ] increases 
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Figure 2. Classification of cells depending on the time constants of peak amplitude decay of Ca 21 spikes. Representative traces of Indo- 
5F signal changes (F/F 0 ; top) and IRIS-1 signal changes (AR/R 0 ; bottom) in cells showing sustained Ca oscillations (S-cell) (A) and damped oscillations 
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(D-cell) (B). Three different concentrations of histamine (1, 3, and 10 uM) were sequentially applied to the same cells with an interval of 20 min. (C) 
Relationships between the histamine concentrations and the inverse time constants of Ca 2+ amplitude decay in S-cells (red) and D-cells (blue). (D) 
Relationships between the histamine concentrations and the Ca 2+ oscillation frequencies observed in S-cells (red) and D-cells (blue). (E-G) 
Comparisons of the Ca 2+ oscillation frequencies between S-cells (red) and D-cells (blue). The histamine concentrations were 1 |jJVI (E), 3 uM (F), and 
10 |j.M (G). (H-J) Comparisons of integrated IRIS-1 signals between S-cells (red) and D-cells (blue). The histamine concentrations were 1 uM (H), 3 uM 
(I), and 10 uM (J). 

doi:1 0.1 371 /journal.pone.008641 0.g002 



through the Ca + influx channels [10]. Because knockdown of any 
single PLC isozyme alone did not significantly decrease the 
amount of the IP3 increase (Fig. 4D), multiple PLC isozymes might 
be activated in response to [Ca j elevation. A significant increase 
in [IP 3 ] was observed in cells treated with PLC[J4KD siRNA#2 
(Fig. 4D), even though the treatment with this siRNA did not 
change the expression levels of other PLC isozymes (Fig. S3). 
Because the other siRNA designed for PLC-p4, PLC[J4KD 
siRNA#l, did not increase [IP3] significandy (Fig. 4D), the 
increase in [IP J may be caused by some side-effects of PLCfMKD 
siRNA#2. PLC-83 knockdown cells using two different siRNAs 
showed significant increases in [IP3] evoked by Ca 2+ increase 
(Fig. 4D). The treatments with these two siRNAs did not change 
the expression level of other PLC isozymes (Fig. S3). PLC-53 may 



Table 1. 


Specific primer sequences for RT-PCR analyses. 








Product size 


Target gene RT-PCR primer sequences 3 (5' to 3 ) 


(bp) 


PLC-pi 


F TGTCTCAGCCCCTTTCTCAC 


111 


R ACAGGAGCACTTGGCGATAC 


PLC-P2 


F GAGGACTTCCCAGAACCTGTC 


116 


R TCCTTCGTCATGTAGGGTTTG 


PLC-P3 


F AAGTCCTTCGACCCCTTCAC 


110 


R TAGATGCCCACCTTCCTGTC 


PLC-P4 


F CCGAACTCGCATGGTTATG 


114 


R TCATACACAGCTATTCTCAAGACAG 


PLC-yI 


F CGACTGGATTCCTCCAACTAC 


108 


R GGCCTGGTTCATCTGCATAG 


PLC-y2 


F TCTGCGCTTTGTGGTTTATG 


109 


R CAGAGGAACGGACCTGAATC 


PLC-81 


F CCAGCAGCACTGAAGCCTAC 


114 


R TGAAAGTATAGCCGTGGTAGATG 


PLC-83 


F GGCCTATGTTAGGGCCTTTG 


110 


R TCTTGGAGGTGAGGGTATGG 


PLC-S4 


F CTGTCGTTTACCACGGACAC 


117 


R TCTCCAGGGACAAGATGACTG 


PLC-8 


F GAAAACACCAGCGATCTTCAG 


107 


R TGCCTCTGGTTGTCCGATAG 


PLC-i; 


F TCGGTGCAGTTATATTCATGTG 


117 


R TCTTCTCTGTGCGTGATAATTC 


PLC-r|1 


F GGAAGGGGAAGTTTCTGATG 


112 


R CCTTATGAAAGATTCCACCTGATG 


PLC-r|2 


F GGACATCATCGAGCAGTTTG 


119 


R GGTGCTCAGGGTTGAAGATG 


P-actin 


F TCGTGCGTGACATTAAGGAG 


110 


R GTCAGGCAGCTCGTAGCTCT 


a F and R denote foreword and reverse primers, respectively. 




doi:1 0.1 371 /journal.pone.008641 0.tOOl 





possess some inhibitory effects on the Ca + -dependent IP 3 
production in HeLa cells. 

Figure 4E (1st component) and 4F (2nd component) show the 
results of PLC isozyme knockdown on the [IP3] increases in 
response to the combination of histamine receptor stimulation and 
Ca 2+ increase. Both the 1st and 2nd components of the [IP3] 
increase after addition of 3 uM histamine plus 2 mil Ca" + were 
impaired by knockdown of PLC-pi or PLC-(34. The increasing 
rate of the IRIS-1 signal of the 1 st component was significandy 
reduced by 48% (PLCpiKD siRNA#l), 38% (PLCpiKD 
siRNA#2), 47% (PLCP4KD siRNA#l), and 51% (PLCP4KD 
siRNA^2) compared with the mean value for the control siRNAs 
mGC and hGC (Fig. 4E). The maximal IRIS-1 signal change of 
the 2nd component was significantly reduced by 33% (PLCpiKD 
siRNA#l), 40% (PLCpiKD siRNA#2), 32% (PLCP4KD 
siRNA#l), and 34% (PLCP4KD siRNA#2) compared with the 
mean value for the control siRNAs mGC and hGC (Fig. 4F). 

A P Y 5 n / / 

1 234121 3 4 e g 1 2 (bp) 




anti-31 anti-(33 anti-B4 anti-y1 anti-53 



P-actin — — — — — — — — — — — — — — 

Figure 3. Identification of PLC isozymes expressed in HeLa cells 
and their specific knockdown using siRNAs. (A) Results of RT-PCR 
analyses using total RNA isolated from HeLa cells. (B) Western blotting 
analyses of endogenous PLC isozymes expressed in HeLa cells using 
isozyme-specific antibodies. Molecular mass markers are shown on the 
right. (C) Western blotting analyses of total lysates prepared from HeLa 
cells treated with PLC isozyme-specific siRNAs. The antibodies used are 
shown at the bottom. The results are representative of at least three 
independent experiments. 
doi:10.1371/journal.pone.0086410.g003 
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Figure 4. Identification of PLC isozymes involved in IP 3 generation in HeLa cells stimulated with histamine. (A) Representative traces of 
IRIS-1 signal changes (AR/R 0 ; top) and lndo-5F signal changes (F/F 0 ; bottom) observed in thapsigargin-treated HeLa cells. The horizontal broken lines 
indicate the baseline levels. (B) Traces of the mean ± SD IRIS-1 signal changes (AR/R 0 ) observed in thapsigargin-treated HeLa cells after addition of 
3 uM histamine (blue circles) or 3 ulVi histamine plus 2 mM Ca 2+ (green circles). (C-F) Effects of PLC isozyme knockdown on the IP 3 increase evoked 
by 3 jiM histamine alone (C), the IP 3 increase evoked by 2 mM Ca 2+ alone (D), the first component of the IP 3 increase evoked by 3 |j,M histamine plus 
2 mM Ca 2+ , and the second component of the IP 3 increase evoked by 3 jiM histamine plus 2 mM Ca 2+ . Data are shown as means ± SD. The numbers 
of cells measured are shown in parentheses. Statistical analyses were performed by one-way ANOVA followed by Scheffe's multiple comparison test. 
*P<0.05, **P<0.01, vs. the values in control siRNA-transfected cells. 
doi:1 0.1 371 /journal, pone.008641 0.g004 



These findings indicate that both PLC-fil and PLC-P4 are 
activated by histamine stimulation accompanied by [Ca *] 
elevation and contribute to the generation of the 1st and 2nd 
components. Delayed inhibition of IP3 production after the 1st 
component of the [IP3] increase was also detected in PLC-Pl and 
PLC-(34 knockdown cells (Fig. S5), suggesting that this inhibitory 
effect was mediated by components upstream of PLC (GPCR 
and/or G proteins) or by both PLC-Pl and PLC-(54 to the same 
degree. 

Roles of PLC-pi and PLC-(34 in the Generation of Ca 2+ 
Oscillations 

We examined how PLC-fil and PLC-fJ4 contribute to the 
generation of Ca 2+ dynamics by measuring [Ca 2+ ] in PLC-fil and 
PLC-(34 knockdown cells stimulated with 3 uM histamine. We 
found that specific knockdown of PLC-fJl resulted in alterations of 
both Ca" + dynamics and IP3 dynamics (Fig. 5A). One of the 
obvious differences was that sustained Ca + oscillations were 
hardly detected in PLC-Pl knockdown cells (Fig. 5B). The mean 



inverse time constants of Ca + spike amplitude decay 
(0.0290±0.0263 s" 1 for siRNA#l and 0.0247±0.01 14 s" 1 for 
siRNA^2) were larger than that in control cells 
(0.0146±0.0093 s -1 , n = 45) and the fraction of cells showing 
damped oscillations was increased (Fig. 5B). In contrast, the mean 
Ca + oscillation frequencies in PLC-Pl knockdown cells were 
similar to that in control cells (37.1 ±8.6 mHz) (Fig. 5C). The 
integrals of the IP 3 signals were significantly reduced in PLC-Pl 
knockdown cells (Fig. 5D). 

When the expression of PLC-P4 was selectively reduced, 
different phenotypes were observed (Fig. 5A). In PLC-P4 
knockdown cells, sustained Ca 2+ oscillations with small inverse 
time constants (less than 10~ 3 s _1 ) were detected in some cells, but 
the overall mean inverse time constants did not differ significantly 
from that in control cells (Fig. 5A and B). The mean frequencies of 
Ca 2+ oscillations in PLC-P4 knockdown cells (21.1 ±7.2 mHz for 
siRNA#l and 17. 0± 10.5 mHz for siRNA#2) were significantly 
smaller than that in control cells (37.1 ±8.6 mHz) (Fig. 5C). The 
integrals of the IP3 signals were also reduced (Fig. 5D), similar to 
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Figure 5. Effects of PLC-pi or PLC-P4 knockdown on Ca 2 and IP 3 dynamics evoked by histamine stimulation. (A) Representative traces 
of lndo-5F signal changes (F/F 0 ; top) and IRIS- 1 signal changes (AR/R 0 ; bottom) in siRNA-treated cells stimulated with 3 uM histamine. The siRNAs 
used are shown on the left. The horizontal broken lines indicate the baseline levels of the IRIS-1 and lndo-5F signals. The vertical broken lines indicate 
the onsets of stimulation. Stacked histograms of the inverse time constants for exponential decay of the Ca 2+ oscillation amplitude (B), Ca 2+ 
oscillation frequencies (C), and integrated IP 3 signals (D) of cells treated with control siRNA (top row), PLC(31KD siRNAs (second row), and PLCp4KD 
siRNAs (third row). The results for the two siRNAs are shown in different colors. The means ± SD are shown at the bottom. The numbers of cells 
measured are shown in parentheses. Statistical analyses were performed by one-way ANOVA followed by Scheffe's multiple comparison test. 
**P<0.01, vs. the values in control siRNA mGC-transfected cells. 
doi:1 0.1 371 /journal.pone.008641 0.g005 



the case for PLC-P 1 knockdown cells. These findings indicate that 
knockdown of PLC-pi increases the time constant of Ca 2+ spike 
amplitude decay, while knockdown of PLC-(34 decreases the Ca 2+ 
oscillation frequency in HeLa cells stimulated with histamine. It is 
important to note that the integrals of the IP3 signals were similarly 
reduced in both PLC-pi and PLC -(54 knockdown cells (Fig. 5D), 
while the temporal patterns of IP3 dynamics were characteristically 
different between these cells (Fig. 5A). In PLC-pi knockdown cells, 
fluctuation of the IP 3 signals was hardly detected during Ca + 
oscillations, while fluctuations of [IP3] synchronized with Ca 2+ 
spikes were obviously detected in PLC-P4 knockdown cells 
(Fig. 5A). These findings suggest that the temporal pattern of IP3 
dynamics, but not the absolute IP 3 concentration, is important for 
determining the pattern of [Ca 2+ ] changes, such as sustained 
oscillations and damped oscillations. 

To confirm the effects of PLC-P 1 or PLC-P4 knockdown on the 
generation of Ca 2+ dynamics, we examined cells transfected with 
PLCP1- or PLCp4-IRES-mRFP constructs (Fig. 6 A) and subjected 
mRFP-positive cells to single-cell imaging of [IP3] and [Ca 2+ ]. 



Overexpression of PLC-pi or PLC-P4 led to different patterns of 
Ca 2+ dynamics (Fig. 6B). PLC-P 1 -overexpressing cells showed 
sustained Ca 2+ oscillations more frequently than control cells 
transfected with the IRES-mRFP construct (Fig. 6B), and the 
mean inverse time constant of Ca 2+ spike amplitude decay 
(0.00349±0.00456 s" 1 , n = 26) was significantly smaller than that 
in control cells (0.0137±0.0100 s" 1 , n= 19) (Fig. 6C). In contrast, 
PLC-P4-overexpressing cells showed Ca 2+ oscillations with higher 
frequency (mean ± SD, 56.8 ± 1 1.6 mHz) than that in control cells 
(mean ± SD, 42.5±4.5 mHz) (Fig. 6D). These findings are 
consistent with those in the RNAi-mediated silencing experiments 
when we make two assumptions as follows. First, histamine- 
induced IP 3 production is mainly operated by PLC-P 1 and PLC- 
P4 in HeLa cells. Therefore, PLC-P4 alone works in PLC-P 1 
knockdown cells and PLC-pi alone works in PLC-P4 knockdown 
cells. No compensation by the other PLC isozymes exists in 
knockdown cells, as supported by the data shown in Fig. S3. 
Second, overexpressed PLC-P isozyme has a dominant activity, 
even though the other type of endogenous PLC-P is expressed in 
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the transfected cells. According to these assumptions, we can 
conclude that PLC-fH activity is involved in determination of the 
time-dependence of Ca 2+ spike amplitude decay, while PLC-P4 
activity is involved in determination of the Ca 2+ oscillation 
frequency in HeLa cells. As previously described in some reports 
[16,17], the overexpression of PLC-|3s alone without G proteins 
did not cause a marked increase in the IP3 production (Fig. 6E). 
Instead, the integral of the IP3 signals was significandy decreased 
in PLC-P4-overexpressing cells to 65% of the value in control cells 
suggesting that an excess amount of PLC-P4 would function 
inhibitory for the IP3 production. The exact mechanism of this 
phenomenon is not clear, but one possibility is that GTP 
hydrolysis by Ga q would be promoted due to an excess amount 
of PLC-f}4 because PLC-Ps are GTPase activating proteins for 
Goc q [12]. When we compared cells showing almost the same 
integral amounts of IP3 production, a phenotypic difference in the 
patterns of [Ca 2+ ] changes was still observed between PLC-P 1- 
overexpressing cells and PLC-P4-overexpressing cells (Fig. 6B). 
Taken together, these findings provide further support for the 
notions that a constant threshold of [IP3] for driving Ca 2+ spikes 
does not exist, and that a temporal pattern of IP3 dynamics is 
critical for the determination of Ca + signaling patterns. 

Discussion 

During the past decade, it has become increasingly evident that 
genetically identical cells can exhibit variability in their cellular 
responses, even in identical environments. Underlying much of 
this variability is the stochasticity in gene expression that can 
produce unique proteomes in individual cells, and this cell 
population heterogeneity is an important factor in processes 
ranging from stem cell differentiation to chemotherapy resistance 
[18,19]. In the present study, we found that individual HeLa cells 
exhibited characteristic and highly reproducible patterns of Ca 2+ 
changes in response to bath applications of histamine, even in the 
same microscopic field of view, as reported in other cell lines 
[5,13-15]. We quantified the patterns of Ca + signals by the Ca 2+ 
oscillation frequencies and the time constants of Ca 2+ spike 
amplitude decay, and found that a difference in histamine 
sensitivity could not account for the cell-to-cell variability in 
Ca 2+ signals. A remarkable finding in this study is that the IP3 
dynamics are characteristically different depending on the patterns 
of Ca 2+ signals (Fig. 1), suggesting that the cell-to-cell variety in 
Ca 2+ signals originates from the upstream of the IP 3 generation 
process in the signal transduction cascade in this cell type. 

Identification of PLC Isozymes Activated in HeLa Cells 
Stimulated with Histamine 

HeLa cells, an epithelial cell line derived from a carcinoma of 
the human uterine cervix, express the histamine H! receptor [20]. 
Ten Ga proteins among the 1 6 mammalian Got proteins and four 
GP proteins among five GP proteins have been detected in HeLa 
cells [21]. RT-PCR analyses positively identified transcripts for all 
12 Gy genes in HeLa cells [21]. In this study, we found that, 
among the 13 PLC isozyme genes, six (PLC-Pl, -P3, -P4, -yl, -83, 
and -s) are expressed in HeLa cells. RNAi-mediated silencing of 
the expressions of the PLC genes showed that at least two isozymes 
in the PLC-P family, PLC-pi and PLC-P4, are actually involved 
in the generation of IP3 production evoked by histamine 
stimulation. Among the three PLC-P isozymes identified in HeLa 
cells, specific knockdown of PLC-P3 did not induce a significant 
decrease in the histamine-induced IP3 production. Recently, 
Adjobo-Hermans et al. showed that PLC-pi and PLC-P4 are 
enriched at the plasma membrane [22], while PLC-P2 and PLC- 



P3 are mainly localized in the cytosol in resting HeLa cells 
following expression of GFP-fused PLC-P constructs. Our findings 
suggest that the intracellular localization of PLC-P isozymes is a 
critical factor for the selection of PLC-P isozymes activated 
downstream of the G protein-coupled histamine H! receptor. 

A Possible Mechanism that Underlies the Generation of 
Cell-specific Patterns of Histamine-induced Ca 2+ Signals 
in HeLa Cells 

The RNAi-mediated silencing experiments and PLC isozyme 
overexpression experiments provide coherent evidence that PLC- 
Pl activity is involved in determination of the time constant of 
Ca + spike amplitude decay during Ca + oscillations, while PLC- 
P4 activity is involved in determination of the Ca 2+ oscillation 
frequency in HeLa cells stimulated with histamine. These findings 
suggest that the relative amounts of PLC-Pl and PLC-P4 have a 
large impact on the pattern of histamine-induced Ca 2+ signals. We 
hypothesize that stochastic variations in the amounts of PLC-P 1 
and PLC-P4 expressed in HeLa cells are involved in the 
generation of the non-genetic cell-to-cell variability in agonist- 
induced Ca 2+ signals (Fig. 7). According to this hypothesis, low- 
frequency sustained Ca 2+ oscillations tended to be observed in cells 
in which PLC-P 1 was dominant, while high-frequency damped 
oscillations tended to be observed in cells in which PLC-P4 was 
dominant (Fig. 7). The former may correspond to the S-cells 
shown in Fig. 2A, while the latter may correspond to the D-cells 
shown in Fig. 2B. Ambler et al. described that sister cells in a 
BC3H-1 clonal cell population displayed unique temporal 
responses in histamine-induced Ca 2+ signals immediately following 
cell division [13], and proposed that the linkage of receptor 
occupancy to Ca 2+ elevation is a functionally unique property for 
individual cells that can be influenced by epigenetic factors. 
Quantification of the individual PLC isozymes and rapid 
manipulation of their expression and/or intracellular localization 
in single HeLa cells will be useful to verify our hypothesis in future 
studies. 

Regulatory Factors for PLC-(3 Activity 

Comparisons of the functional differences between PLC-P 1 and 
PLC-P4 should provide clues toward understanding the mecha- 
nistic basis for the generation of cell-specific patterns of Ca + 
signals. The RNAi-mediated silencing experiments revealed that 
PLC-P 1 only became activated when the histamine stimulation 
was accompanied by [Ca 2+ ] increases in thapsigargin-treated cells, 
whereas PLC-P4 was activated by histamine stimulation with or 
without [Ca "*] elevation. These findings are consistent with 
previous biochemical measurements showing that purified PLC-P 1 
exhibits Ca 2+ -dependent PIP 2 hydrolysis activity with a maximal 
activation at 100-220 nM [Ca 2+ ] [23], while the activity of 
purified PLC-P4 is less sensitive to Ca 2+ and the PLC-P4 activity 
in the absence of Ca 2+ is higher than those of the other PLC-Ps 
[24]. Because feedback regulation by cytosolic Ca + must be 
involved in the Ca 2+ signal formation, the difference in Ca 2+ 
sensitivity for the PIP2 hydrolysis activity between PLC-P 1 and 
PLC-P4 should be critical in determining the pattern of cell- 
specific Ca 2+ signals. 

PLC-P isozymes are regulated by heterotrimeric G proteins, 
whose activities are controlled by the balance of the rate of GTP 
binding, which causes activation, and the rate of GTP hydrolysis, 
which terminates activation [25,26]. Because Ga subunits show 
remarkably slow rates of basal guanine nucleotide exchange and a 
low rate of GTP hydrolysis [27,28], the GTPase catalytic cycle is 
regulated by GPCRs, which accelerate GTP binding, and by 
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Figure 6. Effects of PLC-pl or PLC-P4 overexpression on Ca and IP 3 dynamics evoked by histamine stimulation. (A) Western blotting 
analyses of cell lysates prepared from HeLa cells transfected with PLC-pi -IRES-mRFP (lane 2) or PLC-p4-IRES-mRFP (lane 4). Non-transfected cells were 
used as controls (lanes 1 and 3). (B) Representative traces of lndo-5F signal changes (F/F 0 ; top) and IRIS-1 signal changes (AR/R 0 ; bottom) in 
transfected cells stimulated with 3 ufvl histamine. The plasmid DNAs used to transfect the cells are shown on the left. The horizontal broken lines 
indicate the baseline levels of IRIS-1 and lndo-5F signals. The vertical broken lines indicate the onsets of stimulation. (B-D) Histograms for the inverse 
time constants for exponential decay of the Ca 2+ oscillation amplitude (B), Ca 2+ oscillation frequencies (C), and integrated IP 3 signals (D) in cells 
expressing mRFP (top row), PLC-pl and mRFP (second row), and PLC-p4 and mRFP (third row). The means ± SD are shown at the bottom. The 
numbers of cells measured are shown in parentheses. Statistical analyses were performed by one-way ANOVA followed by Scheffe's multiple 
comparison test. **P<0.01, vs. the values in IRES-mRFP transfected cells. 
doi:1 0.1 371 /journal.pone.008641 0.g006 



GTPase-activating proteins (GAPs), which accelerate hydrolysis. 
PLC-Ps are GAPs for Got q [29] and accelerate the deactivation 
limb of the cycle. In the continued presence of an agonist, it has 
been proposed that activated GPCR stays bound to Ga q through 
multiple cycles of GTP binding and hydrolysis because the rate of 
dissociation of the GPCR from GTP-bound Got q is much slower 
than that of the PLC-P-promoted return of the G protein to the 
GDP-bound state. This phenomenon has been referred to as a 
kinetic scaffold [27,28,30], and the GAP activities of the PLC-P 



isozymes may influence the temporal pattern of the activation of 
each isozyme independently. Recently, it was shown that the 
intracellular localization of PLC-P 1 is regulated by PKC 
phosphorylation [31] and by translin-associated factor X binding 
[32]. In addition to the PIP 2 hydrolysis activity, examination of the 
differences between PLC-P 1 and PLC-P4 in terms of these 
regulatory factors in living cells will be pivotal for understanding 
the molecular basis of their differential roles in the generation of 
Ca z+ signals. 



PLOS ONE | www.plosone.org 



9 



January 2014 | Volume 9 | Issue 1 | e86410 



PLC-fS Isozymes and Cell-Specific Ca + Signals 



Kataoka (Kobe University, Kobe, Japan). Two different prede- 
signed siRNA duplexes targeting each of PLC-(3l, PLC-|33, PLC- 
04, PLC-yl, and PLC-s (Stealth™ RNAi) were purchased from 
Invitrogen. Two Stealth™ RNAi Negative Control Medium GC 
Duplexes (mGCs) and two Stealth™ RNAi Negative Control 
High GC Duplexes (hGCs) (both from Invitrogen) were used as 
negative controls. For knockdown of PLC-53, two original siRNAs 
were synthesized (ST Pharm Co. Ltd.). The sequence of 
PLC83KD siRNA#l (Table 2) was scrambled to minimize the 
sequence homology with any known vertebrate transcripts and 
used to synthesize a negative control siRNA. The nucleotide 
sequences of the siRNAs used in this study are shown in Table 2. 

Plasmid Construction 

The plasmid IRES-mRFP was prepared from IRES2-AcGFPl 
(Clontech Laboratories Inc.) by replacing the AcGFPl cDNA with 
an mRFP cDNA. A full-length mouse PLC-pi cDNA was 
obtained by PCR from mouse PLC P 1 -pBluescriptll KS(+), and 
cloned into the EcoRI and BamHI sites of IRES-mRFP. A full- 
length mouse PLC- (34 cDNA was obtained from mouse PLC (54- 
pBluescriptll KS(+) by Xhol digestion, and inserted into IRES- 
mRFP. 

Cell Culture 

HeLa cells were maintained in DMEM (Nacalai Tesque) 
supplemented with 10% FBS at 37°C under 5% C0 2 . 

RT-PCR 

Total RNA was extracted from HeLa cells using a TRIzol 
Reagent Kit (Invitrogen), and reverse-transcribed using random 
primers and SuperScriptll reverse transcriptase (Invitrogen) 
according to the manufacturer's protocol. cDNA fragments for 
the 13 PLC genes were amplified by PCR using specific primer 
sets (Table 1). The primers were designed using Primer 3 software 
to fullfil the following requirements: (1) the PCR products should 
contain an exon-exon junction to prevent amplification of any 



Table 2. siRNA sequences. 





Name 


Sequences (5' to 3 ) 


GC%" 


Position 1 * 


Control siRNA mGC 


GGUAGGUGAGUGUACAGACGCAAUA 


48 




Control siRNA mGC 


CAGGACGUCUAGCUGUAGUCGCCAU 


56 




PLCpiKD siRNA #1 


GCUACUGGAUCUCAGCCUUGUCAAA 


48 


186 


PLCplKD siRNA #2 


GCCCUCGACCUGAAAUUGAUAACAU 


44 


659 


PLCP3KD siRNA #1 


GGCUUCACUUCGCAUUGCAGCCUUU 


52 


2372 


PLCP3KD siRNA #2 


ACCCGAGACUCAACGAAGUGCUGUA 


52 


766 


PLCP4KD siRNA #1 


GCUGCCAGAUGGUUUCACUGAACUA 


48 


1928 


PLCP4KD siRNA #2 


GCAUGGUUAUGAAUAAUGGACUCAA 


36 


2222 


PLCylKD siRNA #1 


CCCUGCGCUGUAAUGAGUUUGAGAU 


48 


1928 


PLCylKD siRNA #2 


CCUUGUUGACCUCAUCAGCUACUAU 


44 


2199 


PLCeKD siRNA #1 


CGCCACCCUCCAAAGGACUUCAAUA 


52 


1362 


PLCeKD siRNA #2 


CAGCGAUCUUCAGCCUGACCUAGAU 


52 


3792 


Control siRNA scr 


GGCAUAGCCAGAUGGUAGAGAUAAGAA 


44 




PLC83KD siRNA #1 


AGAUG AGCU UCAAGGAG AUCAAGAG AG 


44 


602 


PLC83KD siRNA #2 


UGGUCAACGUGGACAUGAACGACAUAG 


48 


638 



a The GC% column indicates the GC content of each siRNA sequence. 

b The numbers in the position column denote where the sequence is located in the mRNA of the target gene (counted from the first nucleotide of the start codon) of 
human origin. 

doi:1 0.1 371 /joumal.pone.008641 0.t002 
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Figure 7. A possible mechanism that underlies the generation 
of cell-specific patterns of histamine-induced Ca 2 signals. The 

diagrams show the typical IP 3 and Ca 2+ dynamics observed in HeLa cells 
stimulated with histamine. Low-frequency sustained Ca 2+ oscillations 
tend to be observed in cells in which PLC-[S1 is dominant (left), while 
high-frequency damped oscillations tend to be observed in cells in 
which PLC-(54 is dominant (right). 
doi:1 0.1 371 /journal.pone.008641 0.g007 



Materials and Methods 

Materials 

Anti-PLC-Pl, anti-PLC-P3, anti-PLC-P4, and anti-PLC-yl 
rabbit polyclonal antibodies were obtained from Santa Cruz 
Biotechnology Inc. A mouse monoclonal antibody against P-actin 
was purchased from Sigma-Aldrich. A mouse monoclonal 
antibody against PLC-83 [33] was used. Mouse PLC-pi cDNA 
(U85712) and mouse PLC-P4 cDNA (NM_0 13829) were sub- 
cloned and inserted into pBluescriptll KS(+). A rabbit polyclonal 
antibody against PLC-s [34] was a generous gift from Dr. T. 
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contaminating genomic DNA; (2) all known splicing variants 
should be detected; (3) the PCR products should be approximately 
100 bp; and (4) the primers should have almost the same T m 
values to maintain uniform amplification efficiency. Primers for (3- 
actin were designed under similar criteria and used as an internal 
control. The PCR fragments were analyzed by 2% agarose gel 
electrophoresis. 

Western Blotting 

Total HeLa cell lysates were prepared by direct lysis with SDS 
sample buffer (62.5 mil Tris-HCl pH 6.8, 10% glycerol, 2% 
SDS, 5% (3-mercaptoethanol, 0.05% bromophenol blue) and 
subjected to electrophoresis. The separated proteins were trans- 
ferred onto PVDF membranes (Millipore), and the PLC isozymes 
and P-actin were probed with specific primary antibodies. After 
treatment with HRP-conjugated anti-rabbit or anti-mouse sec- 
ondary antibodies (GE Healthcare), the antibody-bound proteins 
were detected with chemiluminescence reagents (ECL: GE 
Healthcare) and a LAS4000mini system (Fujifilm). 

siRNA Transfection 

Fresh culture medium (0.8 ml) was supplied to HeLa cells 
cultured in 35-mm dishes at 30 min prior to transfection. Next, 
1 nl of 20 uM siRNA solution and 1 ul of LipofectAMINE 2000 
(Invitrogen) were independently diluted in 100 ul of Opti-MEM 
(Invitrogen), and the two mixtures were combined and incubated 
for 20 min at room temperature. The mixed solution was added to 
the cells to give a final concentration of 20 nM siRNA. After 
incubation for 8-12 h in a C0 2 chamber, the siRNA solution was 
removed by replacement with 1.5 ml of fresh culture medium. 

Imaging 

For imaging of siRNA-treated cells, HeLa cells at 36-42 h after 
siRNA transfection were transfected with plasmids containing the 
IRIS-1 cDNA using a transfection reagent (TransIT; Mirus). The 
cells were subjected to imaging at 24-30 h after the IRIS-1 
transfection. For imaging of cells expressing recombinant PLC-Ps, 
HeLa cells were transfected with PLCpi-IRES-mRFP, PLCpM- 
IRES-mRFP, or IRES-mRFP plasmids using the TransIT 
transfection reagent at 20-36 h before IRIS-1 transfection. After 
loading the cells with 5 uM Indo-5F AM (Dojindo), imaging was 
performed under a constant flow (2 ml/min) of balanced salt 
solution (20 mM Hepes pH 7.4, 115 mM NaCl, 5.4 mM KC1, 
1 mM MgCl 2 , 2 mM CaCl 2 , 10 mM glucose) as an imaging 
medium at 37°C. For experiments with artificial control of [Ca 
cells were pretreated with 1 nM thapsigargin in the absence of 
extracellular Ca 2+ to deplete the cytosolic Ca 2+ stores by blocking 
sarcoplasmic or endoplasmic reticulum Ca 2+ -ATPase pumps. 
Imaging was performed at 37°C using an inverted microscope 
(IX-71; Olympus) with a cooled charge-coupled device camera 
(ORCA-ER; Hamamatsu Photonics) and a 40 x (NA 1.35) 
objective. A 330-348-nm excitation filter and 460-5 10-nm 
emission filter and a 425-445-nm excitation filter and a pair of 
460-5 10-nm (cyan) and 525-565-nm (yellow) emission filters were 
used for the fluorochromes Indo-5F and IRIS-1, respectively. An 
emission splitter (W-view; Hamamatsu Photonics) was used with a 
fast light source exchanger (DG-4; Sutter Instrument Co.). 
Sequential excitation of IRIS- 1 and Indo-5F was performed using 
a 450-nm dichroic mirror and two excitation filters (425-445-nm 
filter for IRIS-1 and 330-348-nm filter for Indo-5F). Dual 
emission at 460-490 nm (for IRIS-1 and Indo-5F) and 
>520 nm (for IRIS-1) was split with a 460-490-nm filter, a 
long-path 520-nm barrier filter, and two 505-nm dichroic mirrors 
equipped in W-view. Images were acquired at 0.5 or 1 Hz for 



measurements of Ca + oscillations in intact cells and at 0.17, 0.25, 
or 0.5 Hz for measurements in thapsigargin-treated cells. Image 
acquisition was performed with MetaFluor software (Molecular 
Devices). Spectral analyses of Ca 2+ oscillations were performed as 
described previously [35]. 

Supporting Information 

Figure SI Histogram of the time constants for expo- 
nential decay of the Ca 2+ spike amplitude. (A) Ca 2+ spike 
amplitudes were defined as shown in the left scheme. The time 
constants were estimated by fitting the amplitudes with a single 
exponential function, as shown on the right. (B) Histogram of the 
time constants of Ca 2+ spike amplitude decay observed in HeLa 
cells stimulated with 3 |J.M histamine. The arrow indicates the 
threshold for dividing the cells into S-cells and D-cells. 
(PDF) 

Figure S2 Difference in histamine concentration-depen- 
dence of the amplitude of the first Ca 2+ spike and the IP3 
production between S-cells and D-cells. (A) Relationships 
between the histamine concentrations and the peak amplitudes of 
the first increase in Indo-5F signals (F/F () ) in S-cells (red) and D- 
cells (blue). (B) Mean values of the Indo-5F signal changes shown 
in (A). Statistical analyses were performed by one-way ANOVA 
followed by Scheffe's multiple comparison test. *P<0.05, 
**P<0.01. NS: not significant. (C) Relationships between the 
histamine concentration and the integrated IP :i signals observed in 
S-cells (red) and D-cells (blue). Data represent means ± SD. 
Statistical analyses were performed using Student's Hest. 
*P<0.05. 
(PDF) 

Figure S3 PLC isozyme-specific knockdown does not 
affect the expressions of other PLC isozymes. Western 
blotting analyses of total lysates prepared from HeLa cells treated 
with PLC isozyme-specific siRNAs. The isozyme-specific antibod- 
ies used for the western blotting analyses are shown on the left. 
The data are representative of at least two independent 
experiments. 
(PDF) 

Figure S4 Validation of the IRIS-1 signal changes 
observed in thapsigargin-treated HeLa cells. (A and B) 

The maximal IRIS-1 signal changes after addition of 3 uM 
histamine (A) and 2 mM Ca 2+ (B) were plotted against the basal 
CFP intensity of IRIS-1 (top), resting IRIS-1 C/Y ratio (middle), 
and integrated value of Indo-5F change for 10 min (bottom). (C) 
The initial rate of IRIS-1 signal change (1st component) after 
addition of 3 |lM histamine plus 2 mM Ca 2+ was plotted against 
the basal CFP intensity of IRIS-1 (top), resting IRIS-1 C/Y ratio 
(middle), and initial rate of Indo-5F change (bottom). (D) The 
maximal IRIS- 1 signal changes after addition of 3 |iM histamine 
plus 2 mM Ca 2+ (2nd component) were plotted against the basal 
CFP intensity of IRIS-1 (top), resting IRIS-1 C/Y ratio (middle), 
and integrated value of Indo-5F change for 10 min (bottom). 
(PDF) 

Figure S5 [IP3] changes after histamine stimulation 
plus [Ca 2+ ] elevation in thapsigargin-treated PLC-pi and 
PLC-04 knockdown cells. Representative traces of IRIS-1 
signal changes (AR/R 0 ) after addition of 3 uM histamine 
(horizontal open bars) plus 2 mM Ca 2+ (horizontal filled bars) 
observed in thapsigargin-treated HeLa cells transfected with 
control siRNA mGC (left), PLCpiKD siRNAs (middle), or 
PLCP4KD siRNAs (bottom). 
(PDF) 
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